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ABSTRACT 



Or 

<D ] The perturbative QCD fragmentation functions can be applied phenomeno- 

logically as a model for charm and bottom quark fragmentation into heavy- 
light mesons. The predictions by this model on the observables Py and (z) for 
D — D* and B — B* systems are compared with experimental data. 



The dominant production mechanism for mesons and baryons that contain a 
single heavy quark is the fragmentation of the heavy quark, in which light quark- 
antiquark pairs are created out of the vacuum by the color force of the heavy quark and 
then the heavy quark captures the light quarks to form mesons or baryons. However, in 
this process, the creation of light quark-antiquark pairs tells us that the nonperturbative 
q\ ' effects are important and so the fragmentation function, which describes the process, 

cannot be calculated from the first principle or from perturbative QCD. 
Qh! Recently, there are new developments in the fragmentation of heavy quarks into 

heavy-heavy-quark bound states.llm The fragmentation process involves the creation 
of heavy quark-antiquark pair, which implies that the natural scale of the process 
should be of order of the mass of the heavy quark created, and so the process should 
be calculable by perturbative QCD (PQCD). Fragmentation functions for c — > (cc), 
b — > (bb), and b — > (be) in different spin-orbital states have been calculated to leading 
order in a s using PQCD. In addition to the application that these PQCD fragmentation 
functions can predict the production rates of heavy-heavy-quark bound states without 
model dependence, they can also be applied as a phenomenological model to describe 
the charm and bottom quark fragmentation into heavy-light mesons.ffl 

The derivation of the b — ► B c , B* fragmentation functions was given elsewhere.™ 
Here we quote the results: 

rz( 1 — zS^ r 

D- b ^ Bc (z,fi ) = N l (1 _ J [6 - 18(1 - 2r)z + (21 - 74r + 68r V 

-2(1 - r)(6 - 19r + 18r 2 )^ 3 + 3(1 - r) 2 (l - 2r + 2r 2 )z 4 l (1) 



*To appear in proceedings of DPF'94 meeting, Albuquerque, NM (August 1994) 



Fig. 1. (a) The ratio Py = D*/(D + D*) for the D-D* system predicted by the PQCD 
fragmentation functions, and the comparison with the data; (b) the average (z)^ for c — > D* 
and for b —> B* fragmentation versus the scale [i. The experimental measurements from LEP 
Qi = m z /2) and CLEO/ARGUS (ji = 5.3 GeV) are also shown. 



for the So state and 

Tz( 1 — z\ 2 r 

A^ B *(z,/io) = 3iV— -A- i— 2 - 2(3 - 2r)z + 3(3 - 2r + 4r 2 )^ 2 

(1 — (1 — r)zj° 1 

-2(1 - r)(4 - r + 2r 2 )z 3 + (1 - r) 2 (3 - 2r + 2r 2 )z 4 ] (2) 

for the 3 Si state, where N = 2a s (2m c ) 2 | J R(0)| 2 /(8l7rm;?) and r = m c j (rrib + m c ). 

The fragmentation functions in Eqs. ([!]) and (|2]) can also be regarded as two- 
parameter functions with N and r as free parameters. N governs the overall normal- 
ization and r = mi ight /m meson is the mass ratio. We vary the parameter r to study 
the behavior of these fragmentation functions in the limit r = miight/^meson — ► 0. We 
expect that in this limit these fragmentation functions can describe, to certain extent, 
the fragmentation of the heavy quark into heavy-light mesons, e.g., c — > D, D* and 
b — > B, B*. Although in the fragmentation c — ► D, D* and b — * B, B* there are prob- 
ably large nonperturbative and relativistic effects that we have not taken into account, 
the PQCD fragmentation functions with the free parameters N and r can at least 
provide some insights to these systems while precise nonperturbative fragmentation 
functions for charm and bottom are not available yet. 

The ratio Py for the D and D* system is defined as Py = D +* D , , which is a 
measure of the relative population of D* in the production of D and D* mesons. Since 
fragmentation is the dominant production mechanism, the production rates of D and 



D* can be replaced by the corresponding fragmentation probabilities as 



- 5 „ Pc : D : , (3) 



p c ^ D + p c _ 



which is a function of r only. The probabilities P c ~>d an d P c ^d* are obtained by 
integrating D c ^d(z) and D c _>d*(z) over z. The prediction by our fragmentation model 
is shown in Fig. [l|(a). At r = 0, which is the heavy quark mass limit, the ratio Py = 0.75 
is exactly the value given by heavy quark spin symmetry. The average experimental 
value Py = 0.646 ±0.049. For the charm system we choose m c = 1.5 GeV and the light 
constituent quark mass m U)d inside the D and D* mesons to be 0.3 GeV, therefore 
r = 0.167. The experimental data point (r = 0.167, Py = 0.646 ± 0.049) is plotted 
in Fig. |I](a). We found good agreement between the prediction of our fragmentation 
model and the data. Physically, Py < 0.75 means that the production rate of D* is 
less than it should be as given by heavy quark spin symmetry. For the B — B* system, 
with the same value for m Ujd and mj = 4.9 GeV, we have r = 0.058. The value of Py 
at r = 0.058 is about 0.68, which is closer to the heavy-quark-symmetry prediction of 
0.75 than the D — D* system. 

The (z) is the average longitudinal momentum fraction that is transferred from 
the heavy quark to the meson. In terms of fragmentation functions, at a scale \x 
for c — ► D* fragmentation is given by 

/ v M _ Jdz zD c ^ D *(z,fi) 

{ )c ^ D * fdzD c _> D *(z^) ■ W 

The scaling behavior of is given by (z)^ = (z) IMO [a s (iJ,)/a s (iJ,o)}~ 2 ' y ^ b , where 7 = 
— 4Cj?/3, Cf = 4/3, b = (11N C — 2n/)/3, N c = 3, n/ is the number of active flavors 
at the scale fj,, and (z)^ is the value determined at the initial scale fio- Taking the 
inputs: m c = 1.5 GeV, m u ^ d = 0.3 GeV, and [i = m c + 2m u ^ d = 2.1 GeV, we have 

r = 0.167, (z)^ D * = 0.77, and (z)cZd* ^ = 0.50. For b — > B* fragmentation, we choose 
m h = 4.9 GeV, m u>d = 0.3 GeV, and /i = m b + 2m u ^ d = 5.5 GeV, we have r = 0.058, 

{z)t° b * = 0.87, and {z)r~^^ 2 = 0.696. The curves for (z)^ D * and (z)? B „ versus fi 
are shown in Fig. |l|(b). 

Experimental data on (z) are available from the LEP detectors and from CLEO 
and ARGUS detectors. The LEP average {xe) c ^d* = 0.504 ± 0.0133. For the bottom 
quark, only the inclusive hadron production has been measured and {xE)b-*H b = 0.694± 
0.0166. The combined CLEO and ARGUS data is (x e ) c ^d* = 0.648 ±0.043. The scale 
of the measurements is taken to be one half of the center-of-mass energy of the machines. 
These data are shown in Fig. 0(b). Excellent agreement is demonstrated. 

Here, we have demonstrated the applications of the perturbative QCD fragmen- 
tation functions as a phenomenological model to describe the fragmentation of heavy 
quarks into heavy-light mesons and the agreement between the predictions by this 
model and the experimental data. 
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